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We derive asymptotically exact expressions for the thermopower and figure of merit of a quantum 
impurity connecting two noninteracting leads in the linear response regime where the chemical 
potential and temperature differences between the leads are small. Based on sum rules for the 
single-particle impurity spectral function, these expressions become exact at high temperatures 
as well as in the very strongly correlated regime, where the impurity Coulomb repulsion is much 
larger than the temperature. Although modest interactions impede thermoelectric performance, a 
very large Coulomb scale restores the optimal transport properties of noninteracting electrons, albeit 
renormalized to account for the absence of double occupancy in the impurity. As with noninteracting 
electrons, the electronic contribution to the figure of merit is limited only by the spectral broadening 
that arises from the coupling between the impurity and the leads. 
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I. INTRODUCTION 

Nearly twenty years ago, Mahan and Sofo elucidated 
the properties that an ideal thermoelectric material- 
one that can efficiently transform a temperature differ¬ 
ential into a voltage—should have 1 . They found that 
optimal thermoelectric conversion efficiency in the lin¬ 
ear response regime, as determined by a large figure of 
merit, is realized in systems with a transmission function 
T(E) oc 5(E - 6q) exhibiting a Dirac-delta function de¬ 
pendence on energy. In such a system, the electron con¬ 
tribution K e i to the thermal conductance at zero charge 
current vanishes and the figure of merit 


t^el "b Kph 

is only limited by the smallness of K p h■ Here G is the 
electrical conductance, S is the Seebeck coefficient (ther¬ 
mopower), and K e i and n p h are the electrical and phonon 
contributions to the thermal conductance. 

Although realistic systems do not exhibit such trans¬ 
mission, considerable experimental and theoretical work 
has concentrated on exploring thermoelectric properties 
of nearly single-level “impurity” systems such as quan¬ 
tum dots 2-5 or molecular junctions 6-13 . Such systems 
can in principle exhibit a transmission that is strongly 
peaked about this level and hence, promise large figures 
of merit. At the same time, these are small systems 
and their coupling to the leads has a significant effect 
on the transmission function 1~(E), yielding an extrinsic 
broadening T. Nonetheless, the phonon contribution to 
the thermal conductance is expected to be small in these 
systems as compared to that in bulk ones, and there is 
still great interest in exploring their thermoelectric prop¬ 
erties. (Recent studies have stressed the importance of 
the phonon contribution n p h in molecules, however 14,15 .) 

Beyond numerical studies of the thermoelectric perfor¬ 
mance of impurity models based on numerical renormal¬ 
ization group 3,9,16 and the nonequilibrium Green’s func¬ 


tion approach 1 ' ,18 , simple analytic results for thermoelec¬ 
tric coefficients have been found in the “atomic limit”, in 
which the broadening T ->• 0 is taken to be zero, using 
e.g., the sequential-tunnelling approximation for trans¬ 
port 19-21 . In this limit, taking the Coulomb repulsion 
U to be zero, one recovers the idealized, un-broadened 
single-level limit of Mahan and Sofo, and correspond¬ 
ingly, an infinite electronic contribution to the figure of 
merit (i.e., after setting n p h = 0). Turning on interac¬ 
tions, a second satellite level arises due to the interaction 
energy shift U when two electrons are present in the im¬ 
purity. In the atomic limit T -*■ 0 with strong interac¬ 
tions, U » T, Murphy, Mukerjee, and Moore 21 showed 
that the second level becomes unoccupied and hence, ir¬ 
relevant for transport, meaning that the system again 
effectively reduces to a single-level one with a diverg¬ 
ing electronic figure of merit. For intermediate coupling 
strength, U ~ T, both levels are active in transport and 
the figure of merit is not large in general. 

Analogous studies of thermoelectricity in the atomic 
limit of bulk systems have been undertaken in a pair 
of well-known papers by Beni 22 and Chaikin and Beni 23 
(see also Ref. 24). As emphasized by Beni 22 , however, 
it is nontrivial to perturb away from the atomic limit of 
transport since transport formally vanishes when T = 0 
(for them, the limit where the hopping matrix element 
between lattice sites vanishes). For a single-level spin-full 
impurity, the available Hilbert space for transport for¬ 
mally evolves from having three elements, corresponding 
to the three possible occupancies, to an infinite number, 
a large fraction of which must be used in order to have a 
conserving theory of transport 25 . Even when interactions 
are modest, it is challenging to develop a reliable pertur¬ 
bation expansion of transport quantities in powers of the 
coupling matrix elements between the impurity and leads 
(effectively, T). 

In this paper, we use sum rules for the impurity 
electron spectral function to derive expressions for the 
thermoelectric coefficients that are asymptotically exact 
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in two regimes for which the broadening is small, but 
nonzero: First, at high temperatures, greater than 17, T 
and the bandwidth D that characterizes the leads. Sec¬ 
ond, in the very strongly-correlated, “narrow-level” limit 
where U » T » Y and U > D. Our results in this lat¬ 
ter regime reproduce the divergent figure of merit in the 
atomic limit when U -* oo 21 . At the same time, this limit 
is highly singular and we find 

lim ZT(U»T)~ ieo _£ )2 , (2) 

where eo- is the difference between the impurity energy 
and Fermi levels. This result emphasizes the difficulty of 
studying transport in this limit and also the crucial role 
played by the broadening in the thermoelectric perfor¬ 
mance of impurity models with strong interactions. 

We start in Sec. II by introducing the linear-response 
formulae for thermoelectric transport coefficients in 
terms of the single-particle impurity spectral function. 
In Sec. Ill, we introduce two sum rules for this spectral 
function—one that integrates all spectral weight and an¬ 
other which removes the irrelevant weight in the upper 
Hubbard peak—and show how these can be used to de¬ 
rive results for transport that are asymptotically exact 
in the regimes elucidated above. We then calculate ex¬ 
plicit forms for these sum rules for the Anderson impu¬ 
rity model in Sec. IV, and use these results to discuss 
thermoelectric performance in Secs. V and VI. Finally, 
in Sec. VII, we summarize our main results and conclude 
with a discussion of the implications of our results for 
quantum dot and molecular junction systems. 


II. TRANSPORT COEFFICIENTS 


where 

Y = Y l + Y r (6) 

is the sum of the broadenings [assumed in (5) to be purely 
real] at the left and right leads, 

7 e 4r L r fl /r 2 (7) 

is an asymmetry parameter that deviates from unity 

when Tl * Tfl, and 

A(E)^--Y l lmGT(E) (8) 

^ G 

is the spectral function for both spin species a of impu¬ 
rity electrons with retarded Green’s function G ret (E). It 
trivially follows from the analysis in Ref. 27 that (4) is 
also exact with the identification in (5) 9,28 . 

The precise form of the spectral function A(E ) is spec¬ 
ified by the details of the coupling between the impurity 
and the leads. To be specific, suppose 

V = I + h - c -] (9) 

\<LVG 

couples the single-level impurity electrons characterized 
by d(j to the non-interacting v = L,R leads, and 

Hi — (^k if ~ /^i/)hki/cr (19) 

\<LVG 

describes noninteracting electrons with momentum dis¬ 
tribution nk^cr = Ckj^Cki/cr in these leads. For this model, 
A{E) is given by 

A(E) = ZZ( P a+Pb)(b\di\a)(a\d a \b)6(E-E b +E a ), (11) 

g a,b 


Our starting point is the (generalized) Landauer-like 
expressions for the charge 26 

p r°° 

J=— dET(E,T L ,T R )[f L (E)-f R (E)] (3) 

Z7T J -oo 

and heat 

1 r°° 

Jq = — dET(E,T L ,T R )(E-^ L )[f L (E)-f R (E)] 

Z7T J —oo 

(4) 

currents through an impurity connecting two leads. 
fv(E) = {exp[/3j / (if - /j,„)] + 1} _1 is the Fermi function 
at the left [v = L) and right (u = R) leads, with temper¬ 
ature T v e / 3” 1 and chemical potential /i^. Unless speci¬ 
fied otherwise, throughout this paper we set h-k R = 1. 
When the coupling between the impurity and the left and 
right leads are proportional, Meir and Wingreen showed 
that (3) is an exact expression with the identification of 
the generalized transmission function with the spectral 
function A(E) of the impurity electrons 2 ': 

T(E,T L ,T R ) = (n'y/2)T(E)A(E- t i), (5) 


and the broadening is 

r,(U) = 2^^5(F;-ek,)|Vk,| 2 . (12) 

k 


The states |a),|6) in the spectral representation (11) are 
the exact eigenstates of the many-body grand-canonical 
hamiltonian, H\a) = E a \a). Within linear response fiL - 
fJ- R and T l ^ T r , P a = exp (-/3E a )/Z, with Z the grand 
canonical partition function. H includes (9) and (10) 
as well as the as-yet-unspecified impurity hamiltonian. 
We will restrict ourselves in what follows to the linear 
response regime. 

Equations (3), (4), and (5) lead to the following ex¬ 
pressions for the transport coefficients 1 ’ 26 in the linear 
response regime for an energy-independent broadening 
T V (E) = r„ (momentarily restoring h and k R )'- 


G = 


e 2 q r 
16 hk B T 


M 0 , 


(13) 


Kel = 


16 hk B T 2 



Ml) 

M 0 ) 


(14) 
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S = 


and 


ZT = 


1 Mi 
eT Mq ’ 

Ml 


M 0 (M 2 - Ml/M 0 ) ’ 
defined in terms of the integral expression 

M„e dEE n sech 2 (f3E/2) A(E). 

J — oo 


(15) 


(16) 


(17) 


In arriving at these expressions, we have shifted E E+fi 
in the linear-response limit hl - HR - H of the expressions 
(3) and (4) for the currents, meaning that the spectral 
function A(E ) is the one defined in (11). An energy- 
independent broadening is appropriate for leads that ex¬ 
hibit good metallic behaviour with a broad bandwidth, 
generally much larger than T. Had we not made this as¬ 
sumption, the broadening T(E) would have entered the 
integral expression (17). (As we explain below, this situ¬ 
ation can also be dealt with using the methods developed 
in this paper, although the resulting calculations would 
be more complicated.) The above expressions will form 
the basis of our analysis in the remainder of this paper. 

Before closing this section, we briefly comment on the 
major result of Mahan and Sofo who, as noted earlier, 
showed that the electronic contribution to the figure of 
merit would diverge in a system described by a single, un¬ 
broadened energy level. In the language of this section, 
this corresponds to having a Dirac-delta function spectral 
function 


A(E)oc6(E-e 0 ). (18) 


Using this in the above, the thermopower is (restoring 
&b) S = (fc_B/e)(eo/T'). Crucially, using (18) in (17), one 
finds 


M 2 


Ml 

Mo’ 


(19) 


and the electronic contribution to the thermal conduc¬ 
tance vanishes identically at all temperatures , with the 
result that the figure of merit is divergent at all temper¬ 
atures. 

Moving away from the limit (18) of having a single 
un-broadened energy level, a small thermal conductance 
can more generally be understood as expressing a small 
variance in the single-particle energy. Defining 


_ f dE(-)A(E)F(E) 
= / dEA{E)F{E) 


( 20 ) 


which defines the Lorenz number in bulk systems, pro¬ 
vides a direct measure of the variance in the single¬ 
particle energies clustered within T of the Fermi level 
(because of the form of F). It is thus clear why a system 
with a single un-broadened energy level would be opti¬ 
mal, having zero variance in its energy levels; see also 
Eq. (18) in Mahan and Sofo 1 . This result also presages 
our main conclusion: Even though an extremely large U 
can eliminate the energy-level variance due to interac¬ 
tions, unless the broadening T due to coupling with the 
leads can be made smaller than the temperature (not 
currently the situation in experiments 6,29 ), the energy 
variance will be greater that T 2 , and one should only 
expect modest thermoelectric performance. 


III. SUM RULE EXPRESSIONS FOR 
TRANSPORT COEFFICIENTS 

The evaluation of the spectral function A(E) that en¬ 
ters the expressions for the thermoelectric coefficients is 
a challenging many-body problem. In this paper, we pro¬ 
pose an alternative, evaluating instead 

Ml= [°°dEE n A(E ) (23) 


M n (Q c )= / dEE n A(E). (24) 

The first of these is a well-known sum rule, and can 
straightforwardly be evaluated in terms of commutators 
involving electron creation and annihilation operators 
and the hamiltonian (see e.g., Refs. 30,31). Integrating 
the product of E n and (11) leads to (see Appendix A) 

Mj = £({dt,K,M]„}). (25) 


Here [d a ,H] n is a nested commutator with [d a ,H]o = 
d/y the zeroth-order commutator, [d a ,H ]i = [da-,H] the 
first-order commutator, [c? CT ,i7] 2 = \[d a ,H\,H\, and so 
on. 

The second expression, (24), involving an energy cutoff 
H c , is evaluated by projecting out the states with energy 
above H c . Taking T « « U, this amounts to de¬ 

riving an effective operator d a ^ 0 for the lower Hubbard 
peak. The resulting sum rule for the spectral weight con¬ 
tained in the lower Hubbard peak of a single-level impu¬ 
rity model is 


we see that 

M 2 

M 2 --± = M 0 [((E 2 )) F -((E)) 2 F ], (21) 

Mo 

with F = sech 2 (/3if/2). Hence, the ratio 

H = Jp («£ 2 »f - , (22) 


( 26 ) 

<J 

The precise choice of operator d a , f (] depends on the model 
under consideration; its derivation for the Anderson im¬ 
purity model coupled to two leads will be given in the 
next section. We could equally derive a sum rule for the 
spectral weight contained in the upper Hubbard peak, 
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FIG. 1: Schematic plot of the spectral function A(E) at high 
temperatures T » U, T (a) and for strong correlations T < 
T « U (b) relative to the lead Fermi level (denoted here 
by E = 0). Also shown by the dashed lines is the thermal 
weight function sech 2 (/3E/2) that enters the moments M n \ 
the dotted lines show the weight functions—unity and the 
step function 0(f2 c - E )—involved with the sum rules. 

as would be relevant for the situation when the chemical 
potential is close to this peak. In what follows, we will 
assume that the Fermi level is close to the lower Hubbard 
peak, however. 

The sum rule (23) provides a rigorous upper bound on 
M n when n is even: 


show two schematic plots of the impurity spectral func¬ 
tion along with the weight function sech 2 (/3E/2) that 
enters the exact transport moments M n . In Fig. 1(a), 
the temperature is much larger than the Coulomb repul¬ 
sion U as well as T; the majority of the spectral weight 
arises in a Lorentzian centred around the Fermi level of 
width ~ r. A small peak—the weight of which vanishes 
with decreasing UIT —arises at U corresponding to the 
energy of a doubly occupied impurity. In this regime, the 
sech 2 (/3£’/2) function is essentially unity wherever there 
is nonzero spectral weight and one can replace M n by the 
high-temperature asymptotes Mj, with the result that 
the thermoelectric transport coefficients calculated using 
these sum rules become asymptotically exact: 

lim {G, Kei , ZT} [M n \ = {G, /s rf , S, ZT} [M%]. 

T»U,r,D 

(30) 

We formally require that T « D be much larger than 
the lead bandwidth since the spectral function scales as 
A{E) oc E~ 2 for large E, meaning that MJ and 
cannot well-approximate M\ and M 2 unless D <T. 

Fig. 1(b) shows the spectral function in the strongly- 
correlated regime T « T « U. Here two Hubbard peaks 
develop, the lower one centred at some single-particle 
energy scale, eo, and the upper one centred at eo + U. 
When the Fermi level (denoted by the zero in these 
plots) is within ~ T of e.g. the lower Hubbard peak, 
the sech 2 (/3i?/2) function excludes the upper Hubbard 
peak. Because there is negligible spectral weight outside 
the two peaks when U » T 34 , in this regime M n can be 
replaced by M n (H c ), where T « Q c « U and the trans¬ 
port coefficients again admit asymptotically exact sum 
rule values: 


M n <Ml n = 0,2,... VT, (27) 

since A(E) > 0 'iE and sech 2 (x) < 1 Vx. This immedi¬ 
ately allows us to write down upper bounds on the charge 
conductance 

p 2 r 

G < -M 0 t (28) 

8 hT 0 v ' 

and, making use of the fact that M^/Mq > 0 32 , thermal 
conductance K e i at zero charge current: 

*el < 7^ M I- ( 2 9 ) 

Related bounds for the thermal conductance due to 
bosonic excitations (e.g., phonons) were derived in 
Ref. 33. 

Because they involve a ratio of moments, (27) does not 
obviously translate to bounds for the Seebeck coefficient 
S and figure of merit ZT. The usefulness of the sum rules 
is that they provide asymptotically exact approximations 
to the transport coefficients at high-temperatures, in¬ 
cluding the asymptotic limit T » U, T, D, and also in 
the strongly-correlated regime T <T « U. In Fig. 1, we 


lim {G, K e i , S, ZT } [M n \ = 
r«T«u 

{G, K eh S, ZT} [M„(r « n c « U)}. (31) 

Even though (31) is asymptotically exact in the limit T « 
T, it is evident that these results will remain qualitatively 
valid for r < T. 

As noted earlier, we could also consider the thermo¬ 
electric performance when the Fermi level is close to the 
upper Hubbard peak by considering the sum rule for 
this peak. Apart from effects that break particle-hole 
symmetry (enhanced by molecular vibrational modes 18 ), 
the thermoelectric performance will be the same in both 
cases. 

In Figure 2, we plot the regimes in the parameter space 
spanned by U/T, TIT where our results hold. Also shown 
is the regime studied by Murphy, Mukerjee, and Moore 21 , 
as well as that studied using NRG by Costi and Zlatic 3 , 
although NRG can study the entire region shown in this 
figure. Even though it might seem that we are dealing 
with a highly restricted region of parameter space—and 
one that is moreover far from most current experiments 
involving molecular junctions for which T ~ G(10T) 6 ’ 29 — 
we reemphasize that it is only when F « T that a sub¬ 
stantial figure of merit can be realized. We also note that 
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FIG. 2: Plot of the regimes where our results hold (U » T > V 
and T » U,V,D), in comparison to the regimes studied in 
the sequential tunnelling approximation of Ref. 21 (all U, 
T/T -*■ 0) and the numerical renormalization group calcula¬ 
tions of Ref. 3 (U /r = 8, all F). Note that the latter calcu¬ 
lation could have been extended to any region shown above. 
Also indicated are the regimes (r, U « T and F «T,U » T) 
where the electronic contribution to the figure of merit be¬ 
comes large as well as the parameter regime of typical molec¬ 
ular junction experiments. Although the value of U /T is un¬ 
known, F is typically > 10T 6 ’ 29 . 


this regime is well away from the Kondo regime T < Tk, 
since Tk « T 3,35 

In the remainder of this paper, we investigate the im¬ 
plications of (30) and (31) for an archetypal model of 
strongly correlated “impurity” electrons, the Anderson 
impurity model, generalized to include the coupling be¬ 
tween an impurity and two leads. These sum-rule results 
have been derived assuming a constant broadening func¬ 
tion, as is almost always taken to be the case. One can 
also consider energy-dependent broadenings: for simple 
power-law forms such as r(E) = TE l with l and integer, 
the resulting transport coefficients simply involve higher- 
moment sum rules (n -> n + l). For more complicated de¬ 
pendencies, one must resort to the operator-product ex¬ 
pansion technique 36 of Wilson, Kadanoff, and Polyakov. 


IV. SUM RULES FOR THE ANDERSON 
IMPURITY MODEL 

The single-level Anderson impurity model 35 (AIM) for 
an impurity coupled to two leads is 

H = H d + H, + V. (32) 

Here 

Hd = (eo - n) + Uh^fii (33) 

U 

is the impurity hamiltonian and Hi and V are given by 
(10) and (9), the former with = \i, as appropriate in 
the linear response regime. 


A. Sum rules for the entire spectrum 

The results of a straightforward calculation using (25) 
and (32) are 


II 

to 

(34) 

2(e 0 - ^ + Un d /2), 

(35) 


and 

= 2{e 0 - g) 2 + 2{e 0 - g)Und + U 2 nd + 2'Y J \V\ s _ v \ 2 . (36) 

ki/ 

Here n d = T is the total occupancy of the impurity. 
The factor of two in (34) (and elsewhere in the above) 
counts the number of single-particle states: two for the 
single-level AIM with two spins. Taking the V to be 
independent of momentum, 

E|y k ,| 2 = ElK| 2 iVo f D dE 

kis is J D 

= 2ZHVo£|I4| 2 . (37) 

IS 

Here we have taken the lead DOS N a (E) = No per spin to 
be constant (and independent of spin) with D the half¬ 
bandwidth. Making the same assumptions to evaluate 
the broadening (32), it becomes 

r, = 27rV 0 |K| 2 , (38) 

and thus 

Y,\V k „\ 2 = DT/Tr. (39) 

ki/ 

As noted in (21), the combination M - (M^) 2 / Mq = 
Mq[((E 2 ))f = i - ((£J))|, =1 ] provides a measure of the vari¬ 
ance in the single-particle energies, in this case over the 
entirety of the spectral weight since there is no thermal 
weighting function sech 2 ((3E/2). In order to have a small 
electronic contribution to the thermal conductance and 
hence, a substantial figure of merit, the energy variance 
in the system needs to be small. Equations (34)-(36) give 

=v nd( 2 - nd)+ 2 si^i 2 - (4o) 

1V1 0 A ky 

The first term on the right-hand side gives the variance 
in the interaction energy per particle (proportional to 
U 2 times the compressibility), while the second term de¬ 
scribes the variance in the kinetic energy per particle. 
One can understood this result in terms of Fig. 1: with 
two peaks separated by U , the energy variance is clearly 
~ U 2 . The variance in a single band (either one) however, 
is ~ EkJ^I 2 ~ TD. 

At temperatures greater than both T and U, the vari¬ 
ance is of course large, and (36) indicates a correspond¬ 
ingly large n e i. On the other hand, for lower tempera¬ 
tures, T « U , assuming the Fermi level is close to the 
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lower Hubbard peak, the sech 2 (/3£') factor in M n is sen¬ 
sitive only to the spectral weight in that peak. Corre¬ 
spondingly, K e i is only sensitive to the variance ~ TD in 
this peak. This notion can be made precise by projecting 
out the spectral weight contained in the upper Hubbard 
peak from the sum rules, as in (24). We turn to this now. 

B. Projected sum rules for the lower Hubbard peak 

Relevant to the case shown in Fig. 1(b), we now cal¬ 
culate the sum rule for the lower Hubbard peak (LHP), 
obtained by projecting out the doubly-occupied states 
that give rise to the upper Hubbard peak (UHP). Our 
approach here follows closely that used in studies of the 
bulk electron Hubbard model' 51,37,38 and also the AIM 
itself by Schrieffer and Wolff 39 . To remove the spec¬ 
tral weight associated with the UHP, we need a spectral 
decomposition for the impurity electron operator of the 
form 

da — d a , g 0 + d a ,e Q +u • (4-Q 

Here d a ,E describes an annihilation process that lowers 
the impurity energy by ~ E modulo T, i.e., (b\d a ,E |a) = 0 
unless E a - Eb = E + 0(r). Once a suitable expression is 
found for d a eo , one can proceed to calculate the projected 
sum rules given by (26). 

When P = 0, the decomposition is particularly simple 
since the occupancies h a are good quantum numbers, 
and, defining a = -a, 


^<T,eO (l 

(42) 

da,eo+U = 'dad ( j. 

(43) 


When r + 0, we seek a renormalized operator d a such 
that n a is a good quantum number of the bare hamilto- 
nian H, and [h a ,H] = 0. This allows us to undertake a 
similar spectral decomposition as (41) and identify the 
renormalized operator d a ,e 0 for the LHP needed for the 
sum rules. 

Since fi a commutes with H 0 = Ho(d a -*■ d a ,(Aw ->• 
Cki/cr), where 

#o E H d - fi Y, n^ VG (44) 

kv 

is the total hamiltonian less all kinetic energies, [h CT , H] = 
0 implies [Hq,H] = 0. Relating the renormalized opera¬ 
tors to the original operators by a unitary transformation 
(not to be confused with the thermopower S , but in keep¬ 
ing with widely-used notation) 31, 37-40 d a - e s d a e~ s and 
Ckw = e^Ckiyo-e - ' 5 , this last condition amounts to finding 
a transformation operator S that satisfies 

[H, H a ]=0^ [e s He~ s ,H 0 ] = 0. (45) 


Here, and in the remaining, we use a tilde to indicate a 
transformed operator. 

The fact that h a is a good quantum number means 
that—analogous to (41)— we can decompose any oper¬ 
ator O = 0(d a , Cki/cr) in terms of the eigenstates hd = 
Y,cr n<r = U 2, or equivalently, the corresponding energy 
changes E = ±eo,±(e 0 + U) that they effect, modulo cor¬ 
rections o(r): 

O = O ±eo + O±(e 0 +c/) ■ (46) 

0±e operators are related by hermitian conjugacy: 

O-E = & E . (47) 

Equation (41) is a special example of this more general 
spectral decomposition. Because the eigenvalues hd of Hq 
are good quantum numbers, so are the energy changes E 
defined above, with the result that 

[OeM =EO e + 0{T), E = ±e 0 , ±(e 0 + U). (48) 

Equation (48) allows us to determine S and d a eo to a 
specified order in V/ (eo +U), where V is the characteristic 
size of the in-general momentum dependent coupling Vkj, • 
We expand H = exp (S)H exp(-S') as H - H + [5, H] + ••• 
and S = S W + 5( 2 ) + ••• in powers of V/(e 0 + U ) and use 
this in (45). Requiring that S eliminates the iL ± ( eo+ [/) 
terms in the spectral decomposition (46) of H (thereby 
removing the associated spectral weight) and using (48), 
(45) reduces to 

[S^\H d ]=-f eo+u -f. e ^u. (49) 

Here T E are the terms in the spectral decomposition of 
the kinetic energy contribution T = H - Hq to the hamil¬ 
tonian; note that these are the only terms in the hamilto¬ 
nian that can change the impurity energy by ~ ±(eo + U): 
H ± ( eo+l 7) = T ±(eo+u) . Explicitly, 

T eo = XI e ki/n kl /cT + X [ l/ k I Aj {!yfT J fT ( 1 - h a ) + H.c.] (50) 

k vcj k vcT 

and 

T eo+U = X V*„hadtc k „a- (51) 

kva 

Using (51) and (47) in (49), one finds 

S (1) =- l —[f_ eo -u-f to+u \. (52) 

This is essentially 41 just the Schrieffer-Wolff transforma¬ 
tion used by Schrieffer and Wolff to derive the Hondo 
hamiltonian as the effective low-energy description of the 
AIM 39 . Apart from the presence of the additional single¬ 
particle energy scale eo in this expression, it is also the 
same result as obtained for the bulk Hubbard model, with 
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V replaced by inter-site tunnelling; see e.g., Eq. (7) in 
Ref. 40 and Eq. (10) in Ref. 37. 

Having determined the form of we can now ob¬ 
tain an expression for the renormalized LHP operator ap¬ 
pearing in (41). By expanding d a = exp(S)d a exp(-S) = 
d„ + [S (1 >,«k] h— and taking the spectral decomposition 
of both sides of the resulting expression, one finds per¬ 
turbative expressions for the spectral decomposition d a ^ 0 
and d a ^ 0+ u of the original operator; c.f. (41). Using (49) 
and equating terms of the same order in E gives 


da.eg ~ (1 n a )d a + ; ^a^—eg—U ] 

e 0 + U 


-— [fleo-j/, d a -e o-t/] , (53) 

£o + u 


for the LHP operator. Explicit evaluation of this expres¬ 
sion results in 


—(1 n a )d a 


1 


rr (U(r"£u C r) d c 7 I kl , d^ 

eo + U k „ 


TT E + O[V/(e 0 + U )] 2 . (54) 

e o + U kl/ 


In the limit U -* oo, (54) reduces to (42). 

Using (54) to evaluate the first few sum rules (26), 
making use of the fact that there are no doubly-occupied 
states in the transformed basis, d\d\ = 0, one finds after 
some straightforward but laborious algebra 


Mq = 2-nd - l — E [K,<4c k ^)+H.c.]+0(^) 2 , 

£o + U k V CT v eo+u ’ 

(55) 

Mi = (£o - /i)Mo + E + 0(£r) , (56) 

k vcr 


and 


M‘i = (£o - Mi + (2 - nd) E |44i/| 2 

V.v 


ki/cr 


+ 



Vk^’cl,.,da 


■°(£r)- 


(57) 


The analogue of (55) for the bulk Hubbard model has 
been derived in Refs. 31,37,38. Note that e.g. 

E t£,<4*W> = E *WcLa> (58) 

kz/cr kucr 


is purely real in equilibrium, as required to have no charge 
or heat current. Defining 

V = E (59) 

kz/ a 


such that V e V +F is the lead-impurity coupling (9), 
these give the following result for the variance in the 
LHP: 

M 2 

M 2 - =(2 - n d ) E |4U 2 - E - Me*,,*} 

Mq kz/ kizcr 

Comparing the full sum rules (34)-(36) with the pro¬ 
jected ones (55)-(57) and also the variances, (21) and 
(60), one sees that the projection has two primary effects: 
First, it removes spectral weight associated with the UHP 
and the large Coulomb repulsion U. Second, it renormal¬ 
izes the spectral weight from 2 (the total spectral weight) 
to ~ 2-nd: the low-energy spectral weight in the LHP. 
The variance (60) is essentially the variance of a non¬ 
interacting single level , broadened by £ kj/ |U k; ,| 2 and (V); 
the effects of interactions only enter through the spectral 
weight renormalization 2 2-nd and their effect on the 
coherence {d\ c kl/(T ). To the orders in U/(e 0 + t/) shown, 
one can replace the transformed operators d a , c k!y(T in the 
above with the original operators, d a , c kl/(T , mindful, how¬ 
ever, of the restriction that there be no doubly-occupied 
states. 

We now use the above results to give insight into the 
thermoelectric performance of strongly correlated ma¬ 
terials in the high-temperature and strongly-correlated 
regimes. 


V. THERMOELECTRIC PERFORMANCE AT 
HIGH TEMPERATURES T » U 


Using the results from the previous section, we now 
evaluate the exact high-temperature limiting values of 
the thermopower S and figure of merit ZT for the An¬ 
derson impurity model. In the high temperature limit, 
Ud(T » [/,T) -*■ 1. Using this in, (34)-(36), the high- 
temperature values for (15) and (16) become 


S(T -» 00 ) = 
and, further using (39), 

ZT(T ^ 00 ) = 


eo - /r + [//2 Vg 
eT = ~T 


4V J 

v g 


U 2 + 2DT/tt 
Here, we have introduced the gate voltage 


(61) 


(62) 


V g = (e 0 -ti + U/2)/e. (63) 


Taking D = 0 or P = 0, (61) and (62) reproduce the an¬ 
alytic high-temperature limiting values found in Ref. 21 
[specifically, their Eqs. (19) and (21)] using the sequen¬ 
tial tunnelling approximation 19,20 for T ^ 0. From (62), 
one sees that the diverging figure of merit of Mahan and 













Sofo (valid for all temepratures) is realized in the non¬ 
interacting U « 6q and un-broadened limit D,T « eo- 
The high-T limiting form (61) for the Seebeck coefficient 
has the same form as the “Mott-Heikes” approximation 
S(T -* oo) = [fi(T = 0) - fj\/eT —well-known in studies of 
bulk narrow-band systems 23,42 —with the identification 
/i(T = 0) = eo + U /2, which is indeed the zero-temperature 
chemical potential at half-filling, rid = 1. 

Although these high-temperature results are valid in 
a temperature regime well-outside that found in typi¬ 
cal systems, they can serve as useful checks of numerical 
calculations. Any conserving calculation 20 results in a 
transmission function (single-particle Green’s function) 
that satisfies the sum rules and hence, (61) and (62). 
As we have already noted, the sequential tunnelling ap¬ 
proximation used in Ref. 21 satisfies these asymptotic 
results at high temperatures. Our results also give an¬ 
alytic expressions for the high-temperature asymptotes 
found numerically in Ref. 3 using numerical renormaliza¬ 
tion group 43 . 


VI. THERMOELECTRIC PERFORMANCE IN 
THE STRONGLY CORRELATED REGIME U » T. 


In the strong-correlation, narrow-level regime U » T > 
T, (55)-(57) give 


S(U»T>T) = 


1 

eT 


eo+ Al -M+of-eq 

2 -rid \to+u) 


(64) 


for the thermopower and 

(eo - A*) 2 ( 2 - n d ) 2 + 2 ( e o - M)(V t )( 2 - n d ) + (V f ) 2 


ZT=- 


(2 - n d y £ k „ |F k[ ,| 2 + (2 - n d )[(VtV> - <Vt£>] - (Vt>2 
+ ( 65 ) 
for the figure of merit, where we have defined ((5t-/i) 

(^0 E E ~ n)(dlc kv(T ). (66) 

k av 


Note that we are using the bare electron operators in the 
above and (V) is defined using (59). Although (64) and 
(65) are asymptotically exact in the limit where there 
is a large separation between energy scales, T « T « 
U (also the limit where thermoelectric performance is 
optimized), as noted earlier, we expect them to provide 
a reasonable approximation when T <T. 


The thermopower (64) in the strongly-correlated 
regime and intermediate temperatures (T > T) again 
has the Mott-Heikes form, but with a renormalized low- 
energy measure eo + (V^)/(2 - rid) of the single-particle 
energy that does not include the Coulomb energy scale. 
In fact, apart from the renormalization factor 2 - n d ac¬ 
counting for the absence of doubly-occupied states, this 



FIG. 3: Dependence of the Seebeck coefficient (top) and the 
electron contribution to the figure of merit (bottom) on the 
broadening F for different lead bandwidths: D/T = 100 (black 
solid line), D/T = 20 (blue dashed line), and D/T = 10 (red 
dotted line). Curves are shown for U = 100T and eo - [i = T 
although the behaviour shown here is not sensitive to the 
precise choice of values. 


is precisely the single-particle energy of a noninteracting 
electron in an impurity coupled to leads by (9). Consis¬ 
tent with the findings of Ref. 21, by eliminating the pos¬ 
sibility of having a doubly occupied impurity, extremely 
strong correlations have effectively reduced the transport 
problem to that of a single-level. As a result, transport 
is effectively that of noninteracting electrons, albeit with 
a reduced Hilbert space ( n d + 2) 

As with noninteracting electrons, thermoelectric per¬ 
formance in this regime is limited by the broadening 
due to the coupling between the impurity and the lead. 
Specifically, when the lead-impurity coupling is larger 
than the detuning from the Fermi level, (V) > eo — /u, the 
figure of merit (65) is generically of order unity. It is only 
when the coupling becomes much smaller than the detun¬ 
ing, (V) « eo-ft, that one obtains the diverging electronic 
contribution to the figure of merit found in Ref. 21. Out¬ 
side this limit, even though the strong Coulomb repul¬ 
sion has reduced the problem to an effectively single-level 
one, the variance in the single-particle energy (roughly, 
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the greater of ~ Zk ! /I^ki'| 2 and (V^V) - (V 1 ') 2 ) prevents 
the system from attaining optimal thermoelectric perfor¬ 
mance. 

All terms in the denominator of (65) vanish in the 
r ->• 0 or D -*■ 0 limit, meaning that—barring a nonan- 
alytic dependence on either of these variables—the fig¬ 
ure of merit will exhibit the scaling shown in (2) for 
small r ,D. To check this, in Fig. 3 we show (64) and 
(65) as functions of T using Hartree-Fock (HF) theory 
and (39) for a range of bandwidths. Details of the HF 
calculations are given in Appendix B. Although HF is 
not a quantitatively reliable theory when interactions are 
strong, we emphasize that we are only applying HF the¬ 
ory to evaluate the asymptotically exact (conserving) re¬ 
sults (64) and (65), and we expect it to give qualita¬ 
tively reliable results. We are not presenting the results 
of a (non-conserving) HF theory of transport, the re¬ 
sults of which would be completely unreliable. While 
the thermopower is generically of order unity (in units 
of hs/e) for eo - /j, = T, the figure of merit is singular 
at small broadenings and is well-described by the scal¬ 
ing in (2). For comparison with this plot, in molecular 
systems, T ~ (0.01 0.5)eV 6 ’ 29 or, T = (0.4 - 19)T room , 

with T room = 300K. For leads constructed of a good 
metal such as gold, D ~ leV ~ 40T room . 

Fig. 3 emphasizes that either T or the lead bandwidth 
D need to be much smaller than temperature in order to 
reproduce the substantial figure of merit found in Ref. 21 
for U -*• oo. Reinforcing the fact that the thermoelec¬ 
tric transport properties in the very strongly correlated 
regime U » T is equivalent to that of a renormalized non¬ 
interacting system, we note that Fig. 3 closely resembles 
the 17 = 0 Fig. 1 in Ref. 21. 


VII. DISCUSSION 

Much of the interest in studying the thermoelectric 
performance of quantum dots and molecules originates 
in the likely small phonon contribution to the thermal 
conductance in these systems and also the fact that they 
can be well-approximated as single-level systems. For a 
single-level noninteracting impurity system, the figure of 
merit diverges in the limit that the extrinsic broadening 
arising from the coupling to the leads is much smaller 
than the temperature. Turning on a Coulomb repulsion, 
the appearance of a satellite Hubbard peak leads to an 
increase in the variance of the single-particle energy level 
and hence, a diminished figure of merit. 

Somewhat counterintuitively, when the Coulomb re¬ 
pulsion is very strong, however, the new satellite Hub¬ 
bard level becomes unimportant, as it is too high an en¬ 
ergy scale to be relevant for transport. We have used 
asymptotically exact sum rule expressions in this limit 
to show that the thermoelectric transport coefficients as¬ 
sume essentially the same form as for a noninteracting 
system, slightly renormalized to account for the absence 
of doubly occupied states. The thermopower (64) in par¬ 


ticular assumes the well-known Mott- Heikes form, but 
with an effective noninteracting chemical potential. As 
with a noninteracting system, the transport properties 
in this strong interaction regime are limited only by the 
broadening due to the coupling between the leads and 
the impurity. 

This bodes well for using small impurity-type sys¬ 
tems such as quantum dots and molecules to achieve 
substantial thermoelectric efficiencies. Even though the 
bare (un-screened) Coulomb energy scale is very large 
(U » T) in these systems, it is precisely its large size 
that means it does not lead to a diminishment of trans¬ 
port, as would happen if it was of moderate strength, say 
U ~ T . At the same time, the highly singular dependence 
of the figure of merit on the broadening and also the lead 
bandwidth [c.f. ( 2 )] means that it will be essential to 
have a small broadening arising from lead-impurity cou¬ 
pling or to tailor the bandwidth properties of the leads 
in order to realize the potential of impurity systems. 
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Appendix A: Spectral properties of the 
single-particle Green’s function 

In this Appendix, we give an example of the derivation 
of the general result (25) for n = 1: 

/ oo 

dEEA(E) 

oo 

= E E( p a + Pb){a\d G \b){b\di\a){E b - E a ) 

<y a,b 

= [(a\d a H\b)(b\dUa) - (a\Hd a \b)(b\dl\a)} 

o a,b 

+ EE P b [(b\dl\a)(a\lH\b) - (b\di\a)(a\Hd a \b)] 

o a,b 

= E E [{a\d a Hdl\a) - (a\Hd a d ] a \a)] 

o a 

+ EE p o [wiim - (b\diHi\b)] 

o b 

= £ [(d a Hdi) - (HdJl) + (dU*H) - (d ] a Hd a )] 

O 

= £«&&,£]})■ (Ai) 


Appendix B: Hartree—Fock theory for the AIM 
model 

In this Appendix, we describe how to calculate the 
impurity occupation n^, lead-impurity coherence (V) = 


10 


(V*), and (V^), defined in ( 66 ) using Hartree-Fock (HF) 
theory. Within this approximation, (FV) = (V ) 2 and 
hence, these three quantities are sufficient to evaluate 
the figure of merit (65). 

Momentarily generalizing the AIM hamiltonian (32) to 
allow for different chemical potentials and couplings for 
each spin, /x -*• /x^ and V k„ -»■ 14 ^, one has 

n^W = -(^-) ( B1 ) 

V Ofl(j / rp 

and 

(daCywo) = ~(C\i Va d a ) - ~ ’ ( B ^) 

where Cl is the free energy. At the end of our calcula¬ 
tion, we set /Xf = /xj, = /x and likewise with T4w (Bl) is 
a standard thermodynamic identity, while (B2) follows 
from an application of the Hellmann-Feynman formula 
( dil/d\)T,N = (dH/dX) with A = Vk^ a to (32). 

Since it is an effective single-particle theory, the HF 
free energy is given by the trace of the logarithm of the 
HF Green’s function G a (iu) n ) for spin a impurity elec¬ 
trons: Cl = ~/3 -1 Zna eiuin0+ In [-G^(iu n y\, yielding 

« = 4 £ ^'”° + In [ - + eo - Ha + Uria + £ 

P ncr 

(B3) 

Here oj n = (2n + 1 )//3 are Fermi Matsubara frequencies, 
with integer n, and £kc<r = ek v - Ha- The exp(xw rl 0 + ) 
factor with positive infinitesimal 0 + ensures conver¬ 
gence. Making the same assumptions (constant lead 
DOS, momentum-independent coupling Vk^o-) that we 
used to arrive at (38), (B3) reduces to 

fi = -4 £ e iWn0+ ln[- iu n + e 0 ~Ha + Un 9 - (£)hi(f^§)]. 

P ncr 

(B4) 

We will not make the usual approximation in the AIM 
literature of taking D -*■ 00 , in which case ln[(*w n - 
D)/(ioj n + D)] reduces to 7 rsgn(w n ) (see e.g., Sec. 5.2 
in Ref. 35). 

Applying (Bl) and (B2) to (B3) [for the coherence 
(B 2 ) we need to take the derivative of the free energy 
with respect to Vkw before assuming the coupling to be 


constant], one arrives at the results 


= 4e 


D i(jj n 0 


P n iuj n eo + ^a Uria- + ( 4ir ) ln(*“)] +D ) 


(B5) 


and 


{d a ( 


k ucr ‘ 




ioj n 0 


ft n n £kcr) 


V, 


k vcr 


[xcu n - eo + Ha~Un s + (£) ln(^g)] 


(B6) 


Using this last result to evaluate (V) and ( 66 ) and re¬ 
placing the momentum integration with one over energy 
[assuming again, a constant lead DOS as in (37)] gives 


1 


<v> = ~I 


^ n0+ (£) in (m§) 


P n iuj n e 0 + Ha Un a + ( 47 r )l n ( iaj ]] + £)) 


(B7) 


and 




O ioj n 0 


(£)[ 


2D + ioj n In 


( iu n -D \] 
V. iui n + D )\ 


P n iu n -e 0 + [i a - Un s + (£:)ln(^=g) 


, (B 8 ) 


To obtain the results shown in Fig. 3, we solve (B5) 
self-consistently for n a , evaluating the Matsubara fre¬ 
quency sum numerically by summing between n = - 10 4 -*• 
10 4 . To avoid dealing with the convergence factor, we use 

11 ^ 1 

2 P n iu> n -eo + Ha - Un s + ( 4 ^- ) In ( itd " + D ) 

(B9) 

This can be derived from (B5) by re-writing the frequency 
sum as an integral over the complex plane and manipu¬ 
lating the contour of integration; see e.g. Appendix A in 
Ref. 44. In contrast to (B5), for which the convergence 
factor is needed to deal with a summand that vanishes as 
at large frequencies, the summands in both (B7) 
and (B 8 ) vanish as ( iu> n )~ , and we can simply drop the 
convergence factor. The solution of (B9) is used in these 
and we again sum over the Matsubara frequencies, from 
n = - 10 4 ->■ 10 4 , to obtain (V) and (V^). 
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